Devil's staircase for a nonconvex interaction 
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We study ground-state orderings of particles in classical lattice-gas models of adsorption on crystal 
surfaces. In the considered models, the energy of adsorbed particles is a sum of two components, 
each one representing the energy of a one-dimensional lattice gas with two-body interactions in one 
of the two orthogonal lattice directions. This feature reduces the two-dimensional problem to a 
one-dimensional one. The interaction energy in each direction is repulsive and strictly convex only 
from distance 2 on, while its value at distance 1 can be positive or negative, but close to zero. We 
show that if the decay rate of the interactions is fast enough, then particles form 2-particle lattice- 
connected aggregates which are distributed in the same most homogeneous way as particles whose 
interaction is strictly convex everywhere. Moreover, despite the lack of convexity, the density of 
particles versus the chemical potential appears to be a fractal curve known as the complete devil's 
staircase. 



One of the problems of surface physics is to study or- 
derings of atoms (adatoms) adsorbed on crystal surfaces. 
The phenomenon can be modelled by certain classical 
lattice gases |J-§]. Adatoms are assumed to sit on top 
of substrate atoms of the square lattice [QJ^] . Substrate 
deformation causes effective, long-range interactions be- 
tween adatoms. In many cases of interest, the Hamilto- 
nian of such system can be decomposed into two parts; 
one involving two-body interactions between particles lo- 
cated along the x-direction and the other one along the 
perpendicular y-direction of the square lattice. Such 
one-dimensional interactions constitute sometimes cer- 
tain local perturbations of repulsive and strictly convex 
interactions. Knowing the one-dimensional ground-state 
configurations, we can easily construct two-dimensional 
ground-state configurations of the original model, as ob- 
served in |jj and discussed in the end of our paper. There- 
fore, in what follows we deal exclusively with the one- 
dimensional situation. 

In classical lattice-gas models, considered in the sequel, 
every site of a one-dimensional lattice Z can be occupied 
by one particle or be empty. Then, a (infinite-lattice) 
configuration is an assignment of particles to lattice sites. 



We assume that particles interact only through two-body 
forces and to a pair of particles at lattice sites i and j, 
whose distance is r = \i — j\, we assign the translation- 
invariant interaction energy V(r). The Hamiltonian of 
our system in a bounded region A can be then written in 
the following form: 



H A (X) = V ^ 

{{i,j}:{i,j}nA/0} 



■ j\) Si (X) 3j (X), (1) 



where Si{X) assumes value 1 if in the configuration X 
the site i is occupied by a particle and otherwise value 0. 

If we assume that V(r) is positive for every r > 1 and 
strictly convex, that is V(r)+V(r+2) > V(r+1), then we 
obtain a model used by Hubbard to study orderings of 
electrons in quasi-one-dimensional conductors |^||| and 
by Pokrovsky and Uimin to study orderings of monolay- 
ers of atoms adsorbed on crystal surfaces ||. For any 
given particle density p, Hubbard found the energy den- 
sity e(p) of ground-state configurations. He also showed 
that any ground-state configuration has the following re- 
markable property. Let Xi £ Z be a coordinate of the 
i-th particle. Then there exists a sequence of natural 
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numbers dj such that Xi+j — Xi £ {dj, dj + 1} for every 
i £ Z and j £ N, i.e., the distances between the first 
(second, etc.) nearest neighbors can differ by at most 
one lattice constant. Configurations with such property 
arc known as the generalized Wigner lattices or the 
most homogeneous configurations Q|. 

The ground-state phase diagram in the grand- 
canonical ensemble has been calculated heuristically by 
Bak and Bruinsma H while a proof for the related 
Frenkel-Kontorova model has been provided by Aubry 
M and adapted to the lattice-gas model case by Miekisz 
and Radin [0; it is outlined below. 

In the grand-canonical ensemble, with the fixed chem- 
ical potential p, to find the energy density of a ground 
state we have to minimize 

f(p) = e(p) - pp. (2) 

Now, e(p) is differentiable at every irrational p and is 
nondifferentiable at any rational p pJlQ||. However, as a 
convex function, it has a left derivative d~e(p)/dp and 
a right derivative d + e(p)/dp at every p. It follows that 
to have a ground state with an irrational density p of 
particles, one has to fix p(p) = de(p)/dp. For any ra- 
tional p, there is a closed interval of chemical poten- 
tials [d~e(p)/dp,d + e(p)/dp], where the most homoge- 
neous configurations of density p are the ground-state 
configurations. One can show that the sum of lengths 
of these intervals amounts to the length of the interval 
begining at the end of the half-line, where the vacuum 
is the only ground-state configuration, and ending at the 
begining of the hal-line, where the completely filled con- 
figuration is the only ground-state configuration. The 
particle density versus the chemical potential of parti- 
cles, p(p), is constant on each of the above described 
intervals. Moreover, it is a continuous function on the 
real line. The curve p(p) is classified as a fractal one and 
named the complete devil's staircase J8| — |l0|| - 

For any rational p, there is a unique (up to transla- 
tions) periodic ground-state configuration with that den- 
sity of particles. For any irrational p, there are un- 
countably many ground-state configurations which are 
the most homogeneous configurations. It has been shown 
in [jllj that they all belong to one local isomorhism class. 
It means that they cannot be locally distinguished one 
from another. Every local pattern of particles present 
in one ground-state configuration appears in any other 
within a bounded distance, in other words (more for- 
mally), there exists the unique ground-state measure sup- 
ported by them. 

However, it appears that in the models of adsorption 
jl],|]], V(l) is very small (positive, negative or zero) as 
compared with V(2), violating therefore the convexity 
property of interactions. It may seem that the smal- 
ness of V(l) together with the repelling nature of re- 
maining interactions will force particles to form 2-particle 



lattice-connected aggregates, called 2-molecules, allowing 
remaining particles to be farther apart. Then, because 
the interaction between 2-molecules is strictly convex, 
we would obtain a molecular devils'staircase. However, 
as the folowing example shows, one cannot hope to ob- 
tain such a universal result, independent of the details of 
interactions, as in the case of strictly convex interactions. 
Example 1 

Consider the following convex potential V which in some 
units is given by V(l) = 0, V{2) = 4, V{3) = 2, V(A) = 1, 
and V(r) = for r > 5, and two periodic configurations 
(period 17) of particle density p = 8/17 whose elemen- 
tary cells are of the form: [••ooo»»»ooo«»«ooo] 1 and 
[nooitoonootioo o] 2 , where • stands for a particle 
while o for an empty site. The energy per elementary cell 
in the first case is 2V(2) + 3V(4) = 11 while in the second 
case it is 3V(3) + 7V(4) = 13. Thus the configuration 
[.]2 that consists exclusively of 2-molecules looses against 
the configuration that consists of 2-molecules and 3- 
molecules. This remains true if we do not set V(r) = 
for r > 5 but require that V(r) be positive and stricly 
convex from r — 2 on and sufficiently fast decaying from 
r = 5 on. 

Now, in order to get rid of n-molecules, n > 3, we 
restrict the family of considered potentials, demanding 
that V(2) be relatively strong with respect to all other 
V(r)'s. We have arrived at the following result: 

There are no n-molecules with n > 3 in any ground- 
state configuration in a lattice gas ([!]) with an interac- 
tion energy V such that 7(1) > -2V(2) and V(2) > 
(7W + V(l))/2, where W = E£LgV(r). 

The idea of the proof is to start with an arbitrary con- 
figuration containing n-molecules, removing n — 2 par- 
ticles from every n-molecule, rearranging the resulting 
configuration and then putting back removed particles in 
the form of 2-molecules. One can show that if the above 
conditions are satisfied, then the gain in energy caused 
by removing particles is strictly larger than its loss asso- 
ciated with inserting back the particles. 

The following example shows that configurations con- 
taining isolated particles, called atoms, are unlikely to be 
the ground-state ones. 
Example 2 

Consider the following periodic (period 13) configuration: 
[•oo«»oo»oo«« o] 3 . The elementary cell [.]3 con- 
tains a local configuration of two atoms separated by a 
2-molecule. Let us transform this local configuration into 
another one, consisting of two 2-moleculcs only, by mov- 
ing to each other the left atom and the right particle of 
the 2-molecule and the right atom and the right particle 
of the 2-molecule. The resulting elementary cell reads: 
[o • • o o • • o o o • • o] 4 . One can easily calculate the rela- 
tive energy change of the local configurations considered 
above; for a potential described in Example 1 it is neg- 
ative, showing that the original periodic configuration is 
not a ground-state one. 
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We have generalized the above procedure to any lo- 
cal configuration containing two atoms separated by k = 
0, 1, ... 2-molecules. Particles from neighboring pairs are 
moved to each other and k + 1 2-molecules are created. 
Then, it appears that the internal energy of all considered 
particles strictly decreases and the change is bounded 
from above by the interaction energy between two atoms 
taken with the minus sign while the interaction energy of 
considered particles with all the rest does not increase. 
In this way we proved the following statement: 

Consider a lattice gas (|l]) with an interaction energy 
V such that V(l) < 0. Among configurations that do 
not contain n-molecules with n > 3, the lowest-energy 
configurations consist exclusively of 2-molecules. 

The actual proofs of the two results, discussed above, 
are fairly technical and for this reason they are included 
in a separate paper 

Combining the above two results we obtain that 
ground-state configurations of our model contain only 2- 
molecules. It is easy to see that the effective interaction 
between these molecules is strictly convex - it is a sum 
of strictly convex interactions between particles of the 
molecules. The latter observation enables us to formu- 
late a theorem that constitutes our main result. 
Theorem (2-molecule most homogeneous ground-state 
configurations ) 

Consider a lattice gas ([!]) with a nonconvex interac- 
tion energy V and a particle density p < 1/2. If 
-V(2) < V(l) < and V{2) > (7W + V(l))/2, then 
the ground-state configurations are the 2-molecule most 
homogeneous configurations of particle density p. 
Corollary (Molecular devil's staircase) 
In a lattice gas (Q) whose interaction energy V satisfies 
the conditions given in the theorem, the particle density 
versus the chemical potential of particles, p(p-), exibits 
the complete devil's-staircase structure. 

To summarize, we have proven rigorously the exis- 
tence of the molecular complete devil's staircase in one- 
dimensional lattice-gas models with certain nonconvex 
long-range interactions. 

As explained at the very begining, our main result, 
concerned with one-dimensional lattice gases, applies di- 
rectly to a class of two-dimensional lattice gases mod- 
elling the adsorption of adatoms on crystal surfaces. Such 
models have been investigated numerically by Ishimura 
and Yamamoto Q. In particular, they posed the ques- 
tion about the possibility of the presence of a devil's 
staircase in such models. Due to the symmetry of the 
considered Hamiltonians (^), with respect to space rota- 
tions by 7r/2, the two-dimensional ground-state adatom 
configurations are given by the oblique arrangements of 
one-dimensional ground-state configurations determined 
above. Namely, the adatom configurations at lattice lines 
in the x-direction can be chosen as the one-dimensional 
ground-state configurations, with the condition that, go- 
ing in the positive y-direction, the configuration of the 



next lattice line is shifted by one lattice constant in the 
positive rr-direction. Now, knowing the solution of the 
ground-state problem, it would be very interesting to 
investigate the low-temperature stability of such struc- 
tures. 

There is also a numerical evidence [|l7) that certain 
nonconvex interactions appear, as effective two-body 
interactions, in the one-dimensional spinless Falicov- 
Kimball model Jl3|,[l4|]; it is known that this quantum 
model can be transformed into a classical lattice-gas 
model with fairly complicated longe-range and many- 
body interactions 

Another class of systems, where our ground-state 
results are immediately applicable, consists of three- 
dimensional layered systems, like those studied by Fisher 
and collaborators j^J. As in the two-dimensional case, it 
is highly desirable to study the stability of ground-state 
orderings against thermal fluctuations. 
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